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Abstract 



We consider linear stochastic differential-algebraic equations with constant coefficients and 
additive white noise. Due to the nature of this class of equations, the solution must be 
defined as a generalised process (in the sense of Dawson and Fernique) . We provide sufficient 
conditions for the law of the variables of the solution process to be absolutely continuous 
with respect to Lebesgue measure. 
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1 Introduction 

A Differential-Algebraic Equation (DAE) is, essentially, an Ordinary Differential Equation 
(ODE) F(x,x) = that cannot be solved for the derivative x. The name comes from the 
fact that in some cases they can be reduced to a two-part system: A usual differential system 
plus a "nondifferential" one (hence "algebraic" , with some abuse of language) , that is 




x\ = f(xi,x 2 ) 
= g(xi,x 2 ) 



for some partitioning of the vector x into vectors x\ and x%. In general, however, such a splitting 
need not exist. 

In comparison with ODE's, these equations present at least two major difficulties: the first 
lies in the fact that it is not possible to establish general existence and uniqueness results, due to 
their more complicate structure; the second one is that DAE's do not regularise the input (quite 
the contrary), since solving them typically involves differentiation in place of integration. At the 
same time, DAE's are very important objects, arising in many application fields; among them 
we mention the simulation of electrical circuits, the modelling of multibody mechanisms, the 
approximation of singular perturbation problems arising e.g. in fluid dynamics, the discretisation 
of partial differential equations, the analysis of chemical processes, and the problem of protein 
folding. We refer to Rabier and Rheinboldt for a survey of applications. 

The class of DAE's most treated in the literature is, not surprisingly, that of linear equations, 
which have the form 

A(t)x(t) + B(t)x(t) = f(t) , 

with x, f : M + — ► M. n and A, B : M + — > M rixr \ When A and B are constant matrices the equation 
is said to have constant coefficients. Note that these equations cannot in general be split as in 

<dU). 

Recently, there has been some incipient work (Schein and Denk ^2] an d Winkler |14| ) 
on Stochastic Differential- Algebraic Equations (SDAE). In view to incorporate to the model 
a random external perturbation, an additional term is attached to the differential-algebraic 
equation, in the form of an additive noise (white or coloured). The solution will then be a 
stochastic process instead of a single function. 

Since the focus in ^2] and |14j is on numerical solving and the particular applications, some 
interesting theoretical questions have been left aside in these papers. Our long-term purpose is 
to put SDAE into the mainstream of stochastic calculus, developing as far as possible a theory 
similar to that of stochastic differential equations. In this first paper our aim is to investigate 
the solution of linear SDAE with constant coefficients and an additive white noise, that means 

Ax(t) + Bx(t) = f(t) + A£(t) , 

where £ is a white noise and A, B, A are constant matrices of appropriate dimensions. We 
shall first reduce the equation to the so-called Kronecker Canonical Form (KCF), which is easy 
to analyse, and from whose solution one can recover immediately the solution to the original 
problem. Unfortunately, it is not possible to extend this approach to the case of linear SDAE 
with varying coefficients, just as happens in the deterministic case, where several different 
approaches have been proposed. Among these, the most promising in our opinion is that of 
Rabier and Rheinboldt [S]. 

Due to the simple structure of the equations considered here, it is not a hard task to establish 
the existence of a unique solution in the appropriate sense. However, as mentioned before, a 



DAE does not regularise the input f(t) in general. If white noise, or a similarly irregular noise 
is used as input, then the solution process to a SDAE will not be a usual stochastic process, 
defined as a random vector at every time t, but instead a "generalised process", the random 
analogue of a Schwartz generalised function. 

The paper is organised as follows: in the next section we shall provide a short introduction 
to linear DAE's and to generalised processes. In the third section we shall define what we mean 
by a solution to a linear SDAE and in Section |I] we shall provide a sufficient condition for the 
existence of a density of the law of the solution. In the final Section we shall discuss a simple 
example arising in the modelling of electrical circuits. 

Superscripts in parentheses mean order of derivation. The superscript T stands for trans- 
position. All function and vector norms throughout the paper will be L 2 norms, and the scalar 
product will be denoted by (-, •) in both cases. Covariance matrices of random vectors will be 
denoted by Cov(-). The Kronecker delta notation Sij := l{i=f} will be used throughout. 

2 Preliminaries on DAE and generalised processes 

In this section we briefly introduce two topics: the (deterministic) differential-algebraic equa- 
tions and the generalised processes. An exhaustive introduction on the first topic can be found 
in Rabier and Rheinboldt ^Uj, while the basic theory of generalised processes can be found in 
Dawson pQ, Fernique [2], or Chapter 3 in Gel'fand and Vilenkin [3]. 

2.1 Differential- Algebraic Equations 

Consider an implicit autonomous ODE, 

(2.1) F(x,x) = 0, 

where F := F(x,p) : R nxn — > R n is a sufficiently smooth function. If the partial differential 
D p F(x,p) is invertible at every point (xq,pq), one can easily prove that the implicit ODE is 
locally reducible to an explicit ODE. If D p F(xo,po) is not invertible, two cases are possible: 
either the total derivative DF(xq,pq) is onto M n or it is not. In the first case, and assuming that 
the rank of D p F(x,p) is constant in a neighbourhood of (xq,pq), (|2,lj) is called a differential- 
algebraic equation, while in the remaining cases one speaks of an ODE with a singularity at 
(x ,p ). 

A linear DAE is a system of the form 

(2.2) A(t)x + B(t)x = f(t) , t>0, 

where A(t), B(t) G W ixn and f(t) G M. n . The matrix function A(t) is assumed to have a constant 
(non-full) rank for any t in the interval of interest. (Clearly, if A{t) has full rank for all t in an 
interval, then the DAE reduces locally to an ODE.) In the simplest case, when A and B do not 
depend on t, we have a linear DAE with constant coefficients, and an extensive study of these 
problems has been developed. Since we want to allow solutions of DAE in the distributional 
sense, let us precise the definition of a solution. 

Let V be the space of distributions (generalised functions) on some open set [/Cl, that 
is, the dual of the space V = C^°(U) of smooth functions with compact support defined on 
U. An n-dimensional distribution is an element of (T>') n , and, for x = (x±, . . . , x n ) G (D') n 
and <f> G T>, we denote by (x, <fi) = ((zi, 4>), ■ ■ ■ , (x n , (j))) T the action of x on 4>. We will always 
assume, without loss of generality, U =]0,+oo[. 



Definition 2.1 Let f be an n- dimensional distribution on U, and A, B two n x n constant 
matrices. A solution to the linear DAE with constant coefficients 



(2.3) 



Ax + Bx = f 



is an n-dimensional distribution x on U such that, for every test function <p G V, the following 
equality holds: 

A{x,<t>)+B{x,ct>) = {f,ct>) 

The theory of linear DAE starts with the definition of a regular matrix pencil: 

Definition 2.2 Given two matrices A,B£ W nxn , the matrix pencil (A,B) is the function 
A i— > XA + B, for A £ R. It is called a regular matrix pencil if det(XA + B) ^ for some A. 

If the matrices A and B in equation ()2.3)1 form a regular matrix pencil, then a solution exists. 
This is a consequence of the following classical result due to Weierstrass and Kronecker, which 
states that A and B can be simultaneously transformed into a convenient canonical form (see 
e.g. Griepentrog and Marz |I] for the proof). 

Lemma 2.3 Given a regular matrix pencil (A,B), there exist nonsingular n x n matrices P 
and Q and integers < d, q < n, with d + q = n, such that 



PAQ 



Id 
N 



and PBQ 



J 

L 



where Id,Iq are identities of dimensions d andq, N = blockdiag(N\, 
matrix 

fo 1 ... o\ 

1 ••• 



,N r ), with Ni the qi x qi 



Ni 



and J is in Jordan canonical form. 



1 

0/ 



Proposition 2.4 If (A, B) is a regular pencil, x is a solution to the DAE \2. 3\) if and only if 
y = Q~ 1 x is a solution to 

(2.4) PAQij + PBQy = Pf , 

where P and Q are the matrices of Proposition \2.cA 

Proof: The result is obvious, since 1)2. 4|) is obtained from p. 3)1 multiplying from the left by the 
invertible matrix P. □ 
System ()2.4)l is said to be in Kronecker Canonical Form (KCF) and splits into two parts. 
The first one is a linear differential system of dimension d, and the second one is an "algebraic 
system" of dimension q. Denoting by u and v the variables in the first and the second part 
respectively, and by b and c the related partitioning of the vector distribution Pf we can write 
the two systems as follows: 



(2.5) 



















i- 




w 




W 







(2.6) 
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1 + 




1- 




W 




W 




u 



We refer to u as the differential variables and to v as the algebraic variables. 

The differential system has a unique solution once an initial condition, i.e. the value of the 
solution at some suitable test function (fro, is given. The function must have a nonvanishing 
integral (see Schwartz [T^], V.6). It can be assumed without any loss of generality that J °° (fro = 
1. 

On the other hand, system Q2.6JI consists of a number of decoupled blocks, which are easily 
and uniquely solved by backward substitution, without the need of any additional condition. 
For instance, for the first block, 



(2.7) 











( Cl \ 




I- 




1- 




VW 








vw 



a recursive calculation gives the following distributional solution: 



(2i 



<«i,0>=E<pfc^ (fc ~ i) ; 



1. 



,91 



We can thus state the following proposition and corollary: 

Proposition 2.5 Assume (A, B) is a regular matrix pencil. Then, for every vP = (u®, . . . , u9) G 
M d , and every fixed test function (fro with (fro = 1> i/iere exists a unique solution y to \2.J$ 
such that 



(y,</>o)i 



1, 



,«S)e 



Corollary 2.6 Assume (A,B) is a regular matrix pencil. Then, for every u° = (u^, 
M. d , and every fixed test function (fro with J °° (fro = 1, i/iere exists a unique solution x to 12. 3\) 
such that 

(Q~ 1 x,(fr }i = u° i , i = 1, . . . ,d . 



Note that the matrix N is nilpotent, with nilpotency index given by the dimension of its 
largest block. The nilpotency index of N in this canonical form is a characteristic of the matrix 
pencil and we shall call it the index of the equation (|2.3j) . The regularity of the solution depends 
directly on the index of the equation. 



Remark 2.7 Without the hypothesis of regularity of the pencil, a linear DAE may possess an 
infinity of solutions or no solution at all, depending on the right-hand side. This is the case, 
for instance, of 





1 1 

with any fixed initial condition. 



x(t) + 



1 1 




x(t) 



hit) 

Ht) 



2.2 Generalised processes 

As before, let T>' be the space of distributions on an open set U. A random distribution on U, 
defined in the probability space P), is a measurable mapping X : — > (V , B(V)), 

where B(T>') denotes the Borel cr-field, relative to the weak-* topology (equivalently, the strong 
dual topology, see Fernique [2]). Denoting by (X(u>), 4>) the action of the distribution X(u) € V 
on the test function (j) € £>, it holds that the mapping cu i— > (X(u;), 0) is measurable from (O, J 7 ) 
into (R, <B(R.)), hence a real random variable (X,(f)) on (fi,^-", 7-"). The law of X is determined 
by the law of the finite-dimensional vectors ((X, (pi), . . . , (X, (j) n )), 4>% G T>, n G N. 

The sum of random distributions X and 1" on (O, J 7 , T 3 ), defined in the obvious manner, is 
again a random distribution. The product of a real random variable a and a random distribu- 
tion, defined by {aX,4>) := a(X,<j)), is also a random distribution. The derivative of a random 
distribution, defined by (X,(p) := —(X,(j)), is again a random distribution. 

Given a random distribution X, the mapping X: T> — > L°(f2) defined by (j) i— ► (X,(j)) is 
called a generalised stochastic process. This mapping is linear and continuous with the usual 
topologies in X> and in the space of all random variables L (O). Note that we can safely overload 
the meaning of the symbol X. 

The mean functional and the correlation functional of a random distribution are the deter- 
ministic distribution eft i— > E[(X, </>)] and the bilinear form (</>, ■0) i— > E[(X, 0) (X, ?/>)], respectively, 
provided they exist. 

A simple example of random distribution is white noise £, characterised by the fact that (£, (/>) 
is centred Gaussian, with correlation functional E[(£, <?!>)(£, "0)] = 4>(s)ip(s) ds. In particular, 
(£, 0) and (£, ■0) are independent if the supports of <fi and ^ are disjoint. In this paper we will 
use as the base set the open half-line U =]0, +oo[. White noise on U coincides with the Wiener 
integral with respect to a Brownian motion W: Indeed, if (j) is a test function, then 

i>oo 

(2.9) (£,<£)= / 4>{t)dW{t) 

Jo 

in the sense of equality in law. More precisely, the Wiener integral is defined as the extension 
to L 2 (M+) of white noise (see Kuo [7j for a construction of the Wiener integral as extension of 
white noise). Now, integrating by parts in Q2.9|) . we can write 

POO 

^A) = - w(t)j>(t)dt = -(w,j>) , 

Jo 

so that £ is the derivative of the Brownian motion W as random distributions. A random 
distribution is Gaussian if every finite-dimensional projection is a Gaussian random vector. 
This is the case of white noise and Brownian motion. 

Further results on random distributions and generalised stochastic processes can be found for 
instance in the classical papers by Dawson 1 and Fernique |2j. We will also use in Section 3 the 
following facts about deterministic distributions, which apply as well to random distributions. 

The hyperplane TC of T> consisting of those functions whose integral on U is equal to zero 
coincides with the set of test functions which are derivatives of other test functions. Therefore, 
fixing a test function <po € T> such that j v (f>o(t) dt = 1, every cp € T> can be uniquely decomposed 
as (p = A0o + 0; f° r some ip G T> and A = J v <fr(t) dt. 

If / G T>' is a distribution, the equation T = f has infinite solutions (the primitives of /): 
T is completely determined on 7i by {T,ip) = —(f,ip) whereas (T, 0o) can be arbitrarily chosen 
(for more details see Schwartz ^3], II. 4). 



3 The generalised process solution 



Consider the linear stochastic differential- algebraic equation (SDAE) with constant coefficients 

(3.1) Ax + Bx = f + A£ , 

where A and B are n x n real matrices, / is an n-dimensional distribution, A is an n x m 
constant matrix, and £ is an m-dimensional white noise: £ = (£1,... ,£ m ), with £j independent 
one-dimensional white noises. Recall that we will always take U =]0,+oo[ as the base set for 
all distributions. 

Definition 3.1 A solution to the SDAE 

(3.2) Ax + Bx = f + A£ 

is an n-dimensional random distribution x such that, for almost all u € f2, x(u) is a solution 
to the deterministic equation 

(3.3) A±(lo) + Bx(u) = f + A£(w) , 
in the sense of Definition ^. A 

Theorem 3.2 Assume (A, B) is a regular matrix pencil. Then, for every u° = (n®, . . . , vfy £ 
M. d , and every fixed test function c/)q, there exists an almost surely unique random distribution 
x, solution to 4.V. ijj) . such that 

{Q~ l x,(t) Q )i =u° , i = 1, ...,d . 

where Q is the matrix in the reduction to KCF. Furthermore, the solution is measurable with 
respect to the a -field generated by £. 

Proof: For every w £ fi, we have a linear DAE with constant coefficients, given by ()3.3(l . and 
we know from Corollary 12.61 that there exists a unique solution x(oj) S T>' , satisfying 

(Q~ 1 x(uj),(j)o)i = , i = l,...,d. 

The fact that x(uj) is a linear and continuous mapping from T> into R for every w € fi implies 
that to I— > x(w) is measurable from into (P', B(T>')), hence a random distribution (see 

Fernique 0, section III. 4). We still want to prove that the mapping u >— * x(uj) is measurable 
with respect to the u-field generated by the white noise £. To this end, we will explicit the 
solution as much as possible with a variation of constants argument. 

Let P and Q be the invertible matrices of Lemma 12.31 Multiplying (|3.2|) from the left by P 
and setting y = Qx we obtain the SDAE in Kronecker Canonical Form 

(3.4) ({? ° N )y + (i % V -Pf + PK. 

System (|3.4|) splits into a stochastic differential system of dimension d and an "algebraic stochas- 
tic system" of dimension q. Denoting by u and v the variables in the first and the second systems 
respectively, by b and c the related partitioning of the vector distribution Pf, and by S = {o~ij) 



and R = (pij) the corresponding splitting of PA into matrices of dimensions d x m and q x m, 
we can write the two systems as 



(3.5) 





( Ul \ 
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/6\ 


(3.6) 


N 
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+ P 





Fixing a test function (fto with J °° 00 = 1 and a vector u° G M d , we have for the first one 
the distributional stochastic initial value problem 



(3.7) 

where r\ := b + St;. 

Consider the matrix system 

(3.8) 



U + Ju = 7] 

(u,4> ) = u° 



$ + J$ = 
/ °° $(t) • fo(t) * = J d 



whose distributional solution exists and is unique, and it is a C°° matrix function M — > R rfxrf 
(see Schwartz [13], V.6). Define T := ^it. From ipT7j> . it follows that t = Let 

<f>ij(t)(f>(t) = XijMt) +ipij(t) 

be the unique decomposition of the function $>jj ■ (ft £ T> into a multiple of 0o and an element of 
the hyperplane of derivatives TC (see Subsection 12.2)1 . 
Then, 



i=i 



<2 



(3.9) 



Eh<^<M-(E*7^>^)] =E A ^^'^o)-E(%'E^^ 1 ) • 



fc=i 



k=i j=i 



The terms (Tj,<fto) should be defined in order to fulfil the initial condition. Using the decom- 
position 

&i j (t)<fo(t) = 6 ij <fo(t)+$ j (t) 
and applying formula ()3.9j) to (ft = (fto, it is easily found that we must define 



(T J1 (ft ) = u° J + Y,U,Y,^A 1 ) ■ 



k=l 



Therefore, 

d d d d d 

( Ui , <j>)=J2 A u [«? + E (vk, E ^ik 1 )} - E {nk> E ^*7fc 

i=i fc=i £=i fe=i i=i 

d d d d d 

= E x v u °j + E - E A ^ E - E &i*jfc 

j=l fe=l j=l £=1 j=l 

d d d d 

= E ^"S + E E ( E A >^ - ^) 

j=l fe=l £=1 j = l 

Taking into account that 



= / (^(*)0o(«) - SjtMs)) ds , 

JO 

V^W = / {®u(s)<t>{s) ~ ke<t>o(s)) ds , 



we obtain finally 



Li, LL Lb Lb 

E A ^ n i + E E ( / ( E A u^(s)^o(s) - <Ms)<£0)) ds)*^^) 

j=i fc=i £=i ^° i=i 

c( d d pi 

E A ^° + E E ( / ( A<1) ^ - $ <t>){s) d S ) ^(t) 

3=1 k=l t=l J ° 1 

d d „i 

(3.10) = X v u j + E ( / ( A$( ^o - ds • 

7=1 fc=l ^ 



On the other hand, the algebraic part (|3.6f) consists of a number of decoupled blocks, which 
are easily solved by backward substitution. Any given block can be solved independently of 
the others and a recursive calculation gives, e.g. for a first block of dimension q±, the following 
generalised process solution 

qi m 

(3.11) (^,<t3)=E(ca + E^'^ M ) ' J = 1, ■■■<&> </>e£>. 

fc=j £=i 

By (l3~mi) and (|3~TT|) . we have (u,u) = G(£), for some deterministic function G : (2?') m -> (P')". 
Given generalized sequence {rj a } a C (V) m converging to r\ in the product of weak-* topologies, 
it is immediate to see that G(n a ) converges to G(n), again in the product of weak-* topologies. 
This implies that the mapping G is continuous and therefore measurable with respect to the 
Borel ex-fields. Thus, the solution process x is measurable with respect to the cr-field generated 
by f. □ 

Remark 3.3 In the case when 6 = 0, so that the right hand side in \°J. 7| ) is simply S£, it is 
well known that the solution of the differential system is a classical stochastic process which can 
be expressed as a functional of a standard m- dimensional Wiener process. Indeed, we have, in 
the sense of equality in law, from \3.10\) . 



d m 



' 1 '' lib /*0O ft 

(«i^)=EVj+EE (7 (A*0o-*0)(a)<k-* -1 (*)) <r M dWt{t) 

~~i JO Jo J lk 



j=l k=l 1=1 



Fix an initial time to g]0, oo[. Take a sequence {4>o} n C T> converging in T>' to the Dirac delta 
5t , and with supp </>q C [to — ^,to + and let {& n } n be the corresponding sequence of solutions 
to the matrix system hH.H\) . Then, lim^^oo J Q $ n (/>Q = Id ■ lr to oo [(t) a. e. and we get 



d m 



(iH, 0) = E + E E / ( X ^ r ) iki 3 )^ dW e (s 



3=1 



k=l 1=1 " *o 
d m roo 



Uy (ft. f'OQ PS 

EE/ (/ ^-^(w)^'* -1 ^)) 

fe=l £=1 ^° ^° * fc 



where and A are now the limit of their corresponding sequences. 

Now collapsing in the same way (j) to St, which implies the convergence of X to & and that 
of f Q s Qcf) t° Id • l[t,oo[( s ) a - e -> we arrive at 

d d m ~t 

Ut (t) = EM/ + EE / mm-Hs)) ik a M dw e ( S ) . 



k=l l=\ 



Finally, using that the solution to \3. £)) with 5t in place of <f>o is known to be <3?(t) = e J ^ to \ 
we obtain 



no + 



- J ( s ~ to \SdW(s) 



to 



In a similar way we can express the first block of the algebraic part, if c = 0, as 

qi m 



(3.12) 



•±L "° poo 

k=j 1=1 



and analogously for any other block. 



□ 



4 The law of the solution 

In the previous section we have seen that the solution to a linear SDAE with regular pencil 
and additive white noise can be explicitly given as a functional of the input noise. From the 
modelling viewpoint, the law of the solution is the important output of the model. Using the 
explicit form of the solution, one can try to investigate the features of the law in which one 
might be interested. 

To illustrate this point, we shall study the absolute continuity properties of the joint law of 
the vector solution evaluated at a fixed arbitrary test function (j). We will assume throughout 
this section that the base probability space is the canonical space of white noise: O, = T>', 
T = B(T>'), and P is the law of white noise. This will be used in Theorem 14.41 to ensure the 
existence of conditional probabilities (see Dawson pQ, Theorem 2.12). 

Let us start by considering separately the solutions to the decoupled equations ()3.5|) and 
(J3.6|) , From the explicit calculation in the previous section (equation (|3.1U|) for the differential 
part and equation (|3.11|) for the first algebraic block), we get that for any given test function <f> 
the random vectors (n, <fi) and (v, <ft) have a Gaussian distribution with expectations 

d d „t 

E[(ui,4>)] =E A ^i+E(^'( / (A*0o - *0)(«) d« (<))..) > i = l,-..,d, 
j=i k=i Jo 1 

E[(v i ,4>)]=Y,(ckA {k - i) ) , i = l,-..,qi, 

k=i 



and covariances 
(4 



.1) Cov((u,</>» = V / [V / (A*0o-^)W<fe-* _1 

x [ V / (A$0 O - $<f>)(s) ds ■ $- l (t)) jk a M 
k=i Jo 



{t))ik&kl 



ds 



l,...,d 



and 
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(4.2) 



xcov((e,0),(e,0),...,(e,0 {9i " 1} )) 



/ Pi 92 
92 93 

W 



Pqi-l 9qi\ 
Pit 







0/ 



where denotes the i-th row of the matrix R and Cov((£, #),..., (£, is a square matrix 

of dimension mq\. We refer the reader to |Hj for a comprehensive study of multidimensional 
Gaussian laws. 

For the differential variables u alone, we are faced with a usual linear stochastic differential 
equation (see Remark 13. 3|) , and there are well-known results on sufficient conditions for its 
absolute continuity, involving the matrices S and J (see e.g. Nualart |S]). 

For the algebraic variables v, their absolute continuity depends in part on the invertibility 
of the covariance matrix of the white noise and its derivatives that appear in (|4.2|) . We will 
use the following auxiliary result concerning the joint distribution of a one-dimensional white 
noise and its first k derivatives. This is a vector distribution with a centred Gaussian law and 
a covariance that can be expressed in full generality as (cf. Subsection 2.2) 

cov((c,<A),...,(e (fc) ,<A))^ = RE r ' • — 



(4.3) 



( _i)^"ll0((*+;)/2)| 



where Re means the real part. We can prove the absolute continuity of this vector for k < 3. 

Lemma 4.1 For all (ft € V— {0}, and a one- dimensional white noise £, the vector ((£, £, £, £ ), (ft] 
is absolutely continuous. 



Proof: The covariance matrix of the vector ((£,£,£, £,),(ft) is 







\ 

" 2 







V o 



whose determinant is equal to 



det 



• det 



Both factors are strictly positive, in view of the chain of strict inequalities 



(4.4) 



These follow from integration by parts and Cauchy-Schwarz inequality, e.g. 



< 



and the inequality is strict unless <ft = K<p for some K, which implies = 0. □ 

The proof above does not longer work for higher order derivatives and we do not know if the 
result is true or false. 

Consider, as in the previous section, only the first algebraic block, and assume momentarily 
that its dimension is q\ = 2. From (|4.2[l . the covariance matrix of the random vector {(v i, V2), 4>) 
is 



2 IMI 2 + 



2 IU II 2 

\\P2\\ 



with determinant 



|P1|| 2 ||P2|| 2 



2 (Pl,P2) 

(Pl,P2> 2 ) + 




2 \\n II 4 
l|P2|| 



Hence, assuming (p ^ 0, we see that the joint law of (v±, 4>) and (v%, 4>) is absolutely continuous 
with respect to Lebesgue measure in M 2 if p2 is not the zero vector. When \\p2W = but \\pi\\ 7^ 0, 
then (v2,4>) is degenerate and (i>i,0) is absolutely continuous, whereas 1 1 yc»2 1 1 = ||pi|| = makes 
the joint law degenerate to a point. 

This sort of elementary analysis, with validity for any test function 0, can be carried out for 
algebraic blocks of any nilpotency index, as it is proved in the next proposition. Let us denote 
by £(k) the subset of test functions <f> such that the covariance Cov((£, 4>), . . . , (^ k ~ 1 \ </>)) is 
nonsingular. With an m-dimensional white noise, the covariance is a matrix with (fe + 1) 2 square 



m x m blocks, where the block (i,j) is Re 



-1) 2 



»+i)/2)||2 



times the identity /„ 



Proposition 4.2 Let (v\, . . . ,v qi ) be the generalised process solution to the first block of the 
algebraic system \2. 61) and r the greatest row index such that \\p r \\ 7^ ; and fix 4> £ 
Then ((i>i, • • • , v r ), 4>) is a Gaussian absolutely continuous random vector and ((u r +i, . . . , v qi ), 4>) 
degenerates to a point. 



Proof: We can assume that c = 0, since the terms Ylk=j( c k> 
additive constants. Then we can write 



({viA) 


y 


\{Vqx,4>) 





P2 



P2 
PS 



Pqi-l 
Pqi 



\P 



qi 



Pqi\ 




0/ 



J ') in (|2.6() only contribute as 



If r is the greatest row index with ||p r [| 7^ 0, it is clear that the qi x mq\ matrix 



fPi 

P2 



\P 



91 



P2 
P3 





Pqi-1 
Pqi 





Pqi\ 




0/ 



has rank r. The linear transformation given by this matrix is onto W x {0} qi ~ r . From this 
fact and the absolute continuity of the vector ((£, (j>), . . . , (^ qi ~ 1 \ 4>)), it is immediate that the 
vector ((vi,(p), . . . , (v r ,cp)) is absolutely continuous, while ((v r +i, </>),...,(«, 
a point. 



qi ,4>)) degenerates to 
□ 



Let us now consider the solution x to the whole SDAE (|H.2j) . We will state a sufficient 
condition for the absolute continuity of (x,(ft), (ft G T>. The following standard result in linear 
algebra will be used (see e.g. Horn and Johnson jS], page 21). 

Lemma 4.3 Let the real matrix M read blockwise 

'A B 



M 'C D 

where A G R dxd , B G R dxq , C G R qxd , D G R qxq and D is invertible. Then the d x d matrix 

A-BD^C 

is called the Schur complement of D in M and it holds that 

det M = det D ■ det{A - BD^C) . 

□ 

A natural application of this lemma is in solving a system of linear equations: 

Ax + By = u | 
Cx + Dy = vj 

where x, u G R d , y, v G R q . We have 

(4.5) u = (A - BD~ l C)x + BD~ l v 

and, if M is in addition invertible, the solution to the linear equation is given by 

x = (A - BD~ x C)~ x (u - BD~ l v) 

y = D~ l {v - C{A - BD- l C)- l (u - BD~ 1 )) . 

We now state and prove the main result of this section. 

Theorem 4.4 Assume (A, B) is a regular matrix pencil and that the matrix A of equation h'J.l\) 
has full rank, and call r the nilpotency index of the SDAE \°J.1\) . Then the law of the unique 
solution to the SDAE \3.1)) at any test function (ft G £(r) is absolutely continuous with respect 
to Lebesgue measure on R a . 

Corollary 4.5 Under the assumptions of Theorem \4-4[ if the nilpotency index is r < 4, then 
the law is absolutely continuous at every test function (ft G T> — {0}. 

Proof of Theorem \4-4\ ft w ^ be- enough to prove that the random vector ((u,v),(ft), solution 
to (|3.5|) and (|3.6|) . admits an absolutely continuous law at any test function (ft E £, since the 
solution to the original system is then obtained through the non-singular transformation Q. 

We shall proceed in two steps: first we shall prove that (v, (ft) admits an absolutely continuous 
law, and then that the conditional law of (u, (ft), given (v, (ft), is also absolutely continuous, almost 
surely with respect to the law of (v, (ft). 

Step 1: We can assume c = in (J3.6|) . By Proposition 14.21 the solution to any algebraic 
block is separately absolutely continuous. Assume now that there are exactly two blocks of 
dimensions q\ and q2, with q<i < q\; the case with an arbitrary number of blocks does not pose 
additional difficulties. 



As in Proposition 14.21 we have 



((vi,<f>), {v qi ,(f>}, {v qi+1 ,(f>), (v qi+q2 ,(f>}) 

( Pi P2 

P2 P3 



(4.6) 



Pqi-l Pqi\ 

Pgi 



Pqi 













Pqi+1 Pqi+2 Pqi+q2 

Pgi+2 Pqi+3 ' ' ' Pgi+<J 2 











<3l+<?2 















0/ 



Since the {q\ + ga) x m matrix R = (pi, ... , p gi + g2 ) has, by the hypothesis on A, full rank 
equal to q\ + q2, the transformation defined by (|4.6j) is onto R q . From the absolute continuity 
of the vector ((£, 0), . . . , (f fa" 1 ), 0)), we deduce that of (u, d>). 

S'tep 2: Since the n x m matrix A has full rank we can assume, reordering columns if 
necessary, that the submatrix made of its first n rows and columns is invertible, and that 



A 



ABE 
CDF 



where A G R dxd , B G R dxq , C G R qxd , D G R qxq , E G R<"< ("»-»), F G R9* ("»-"), with invertible 
D. Let us define 

W := (Vi, ... ,V q ,Vi+V2,.. ■ ,V qi -i + V qi ,V qi+ i + V qi+ 2, ■ ■ ■ , Vq 1+q2 -l + %+ ?2 , £n+l , ■ ■ ■ , Cm) • 

We can write then 

G2 



(4.7) 



G 3 



where G\ is the matrix in (|4.6|) . 

/ Pi 



G 2 



Ga '.- 



(e d+q+l ••• N 



\ e m 



i ^5 



/ P91+1 ••• 

\Pgi+? 2 -l " ' 

/ei ••• 0\ 

\ed ••• 0/ 



and (ej)j = 

By the invertibility of D and the fact that the rows p qi and p qi + q2 have at least one element 
different from zero, it is easy to see that the matrix in (|4.7j) has itself full rank and, therefore, 
that the vector ((u>, £1, . . . , ^), </>} has an absolutely continuous Gaussian law. 



Using (|4.5jl . we obtain 
(4.8) 



/ v 1 +v 2 \ 

V2 + * 3 







( Ul \ 




+ J 


: = (A - BD^C) 


\Ud) 




\u d ) 



+ BD 



-i 



J Qi 



+ (E-BD^F) 




Vq 1+ 1 + Vq 1+ 2 
V qi +2 + Vq 1+ 3 

V V 11+12 J 

It is obvious that both Definition 13,11 and Theorem 13.21 continue to hold true with any 
generalised process 9 in place of the white noise £ in the right-hand side. From Theorem 13 . 21 we 
have in particular that the solution u to the differential system (|4.8j) is a measurable function 
G: (V') m -» (V') d of its right-hand side 9, that means, it = G(0). Let 

p : B{V) xV — ► [0, 1] and q : B(V) xV — ► [0, 1] 

be conditional laws of u given w, and of 9 given w, respectively. That means that 

P({u £ B}n{w £ C}) = I p(B,w)(j,(dw) 

Jc 



and 



P({9 G B} D {w G C}) = / q(B,w)n(dw) , 

Jc 



for any B,C £ B(T>'), where \x is the law of w. 

For every w £ V, let Z w be a random distribution Z w : f2 — > D' with law P{Z W £ B} 
q(B,w). 

Then 

/ P{G(Z W ) £ B} fi(dw) = [ P{Z w £G- 1 (B)}v(dw) 
Jc Jc 

= [ q(G- l (B),w)fi(dw) = P({9 £G- 1 (B)}n{w £ C}) 
Jc 

= P{{G(9) £ B} n {w £ C}) = P({u £ B} n {w £ C}) 
p(B, w) fi(dw) 



c 

Therefore P{G(Z W ) £ B} = p(B,w) almost surely with respect to the law of w, for all 
B £ B(T>'). We have proved that if the right-hand side of the differential system has the law 
of 9 conditioned to w, then its solution has the law of u conditioned to w. It remains to show 
that this conditional law is absolutely continuous, almost surely with respect to the law of w. 

Now, for each w, we can take Z w as 

Z W = (A- BD- l C)i] w + a w 

where a w is a constant ci-dimensional distribution, and (rj, 4>) is, for each (ft £ V, a Gaussian el- 
dimensional vector. This random vector is absolutely continuous: Indeed, its law is that of the d 
first components of the m-dimensional white noise (£i, • • • , £ m ) conditioned to lie in an (m — d)- 

1/2 

dimensional linear submanifold. Let L Wj< p be its covariance matrix. Then (r],(j)) = LJ,((,cp), 
for some d-dimensional white noise ( = (£i, • • • , Cd)- 



Consider now the (ordinary) stochastic differential equation 



a w + (A- BD^C)L 




(4.9) 



By hypothesis, the Schur complement A — BD~ 1 C is non-singular, and therefore the matrix 
(A — BD~^C)L l J^. is itself non-singular. But in the situation of (|4.9|) . it is well-know that 
the solution ((n 1; . . . ,u d ), </>) is a stochastic process with absolutely continuous law for any test 
function (ft ^ 0- 

We conclude that the law of (u,<f>) conditioned to (w, ((>), which coincides with the law of 
(u, eft), is absolutely continuous almost surely with respect to the law of w. This is sufficient to 
conclude that ((ui, . . . , Ud, v\, . . . , v q ), (ft) has an absolutely continuous law, which completes the 
proof. □ 



5 Example: An electrical circuit 

In this last section we shall present an example of linear SDAE's arising from a problem of 
electrical circuit simulation. 

An electrical circuit is a set of devices connected by wires. Each device has two or more 
connection ports. A wire connects two devices at specific ports. Between any two ports of a 
device there is a flow (current) and a tension (voltage drop). Flow and tension are supposed to 
be the same at both ends of a wire; thus wires are just physical media for putting together two 
ports and they play no other role. 

The circuit topology can be conveniently represented by a network, i.e. a set of nodes and 
a set of directed arcs between nodes, in the following way: Each port is a node (taking into 
account that two ports connected by a wire collapse to the same node), and any two ports of 
a device are joined by an arc. Therefore, flow and tension will be two quantities circulating 
through the arcs of the network. 

It is well known that a network can be univocally described by an incidence matrix A = (a^ ). 
If we have n nodes and m arcs, A is the m x n matrix defined by 

{+1, if arc j has node i as origin 
— 1, if arc j has node i as destiny 
0, in any other case. 

At every node i, a quantity di (positive, negative or null) of flow may be supplied from the 
outside. This quantity, added to the total flow through the arcs leaving the node, must equal 
the total flow arriving to the node. This conservation law leads to the system of equations 
Ax = d, where Xj, j = 1, . . . , n, is the flow through arc j. 

A second conservation law relates to tensions and the cycles formed by the flows. A cycle is 
a set of arcs carrying nonzero flow when all external supplies are set to zero. The cycle space 
is thus kerA C W 1 . Let B be a matrix whose columns form a basis of the cycle space, and let 
c £ W 1 be the vector of externally supplied tensions to the cycles of the chosen basis. Then we 
must impose the equalities B T u = c, where Uj, j = 1, . . . , n, is the tension through arc j. 

Once we have the topology described by a network, we can put into play the last element 
of the circuit modelling. Every device has a specific behaviour, which is described by an equa- 
tion <p(x,u,x,u) = involving in general flows, tensions, and their derivatives. The system 



<&(x, u, x, u) = consisting of all of these equations is called the network characteristic. For in- 
stance, typical simple two-port (linear) devices are the resistor, the inductor and the capacitor, 
whose characteristic (noiseless) equations, which involve only their own arc j, are ua = Rxj, 
Uj = L±j, and xj = Ciij, respectively, for some constants R,L,C. Also, the current source (xj 
constant) and the voltage source (uj constant) are common devices. 

Solving an electrical circuit therefore means finding the currents x and voltage drops u 
determined by the system 

Ax = d 
B T u = c 
$(x, u, x, ii) = 

Example 5.1 Let us write down the equations corresponding to the circuit called LL-cutset 
(see \1 j) /. pag. 60), formed by two inductors and one resistor, which we assume submitted to 
random perturbations, independently for each device. This situation can be modelled, following 
the standard procedure described above, by the stochastic system 



(5.1) 



x\ = -x 2 = x 3 

U± - U 2 + U 3 = 

ui = Li±i + ri£i 

U2 = L 2 X 2 + T2^2 

u 3 = Rx 3 + t 3 £ 3 



where £i, £ 2 , £3 are independent white noises, and t\,t 2 , t 3 are non-zero constants which measure 

the magnitude of the perturbations. With a slight obvious simplification, we obtain from \5.1\) 
the following linear SDAE: 
(5.2) 

/0 \ /nA / R- 1 -R- 1 1 0\ /uA / -^R-^ 

ii 2 -RT 1 R- 1 1 u 2 _ ^R' 1 

L x x x + -1 xt -n 

\0 L 2 ) \x 2 ) \ -10 0/ \xj \ -r 2 / 



Let us now reduce the equation to KCF. To simplify the exposition, we shall fix to 1 the values of 
Ti, R and Lj. (A physically meaningful magnitude for R and Li would be of order 10~ 6 for the 
first and of order 10 4 for the latter. Nevertheless the structure of the problem does not change 
with different constants.) The matrices P and Q, providing the desired reduction (see Lemma 
HO)) , are 



f 



P 



V-i 



1 


0/ 



Q 



V-i 



Indeed, multiplying by P from the left and defining y = Q 



_3 

4 
_i 

4 

1 

2 
1 
2 

L .r, 



(5.3) 



(I 0\ 

10 



\0 0/ 



/| o\ 
0100 
0010 

\o 1/ 



yit) 



f-n 

-Tl 



V 



T 2 
-T2 






-1\ 




0/ 

we arrive to 
-r 3 \ 

~T 3 


T3 ) 



We see that the matrix N of Section has here two blocks: A single zero in the last position 
(y±) and a 2-nilpotent block affecting y 2 and y 3 . We have therefore an index 2 SDAE. From 
Propositions \4-2\ and Theorem \4-4\ we can already say that, when applied to any test function 



(f) ^ 0, the variables 2/4, y 2 and y\, as well as the vectors (2/1,2/2) and (2/1,1/4), will be absolutely 
continuous, whereas 2/3 degenerates to a point. 

In fact, in this case, we can of course solve completely the system: The differential part is 
the one- dimensional classical SDE 

(5.4) yi + -2/1 = -ti£i + t 2 £ 2 - T36 , 
and the algebraic part reads simply 

2/3 + 2/2 = -Tl6 - T26 - T36 

(5.5) { 2/3 = 

2/4 = T3C3 • 

XYie solution to \5.'J\) can thus be written as 

2/1 (i) = e"^')/ 2 L(t ) + / e-( s - i0 )/ 2 (-r 1( iWi + r 2( i^ 2 - r 3 dW3)(s) 

2/2 = -nil - t 2 & - T36 
2/3 = 
2/4 = T36 , 

where Wi, W2, W3 are independent Wiener processes whose generalised derivatives are £1, £ 2 
and £3. Multiplying by the matrix Q we finally obtain the value of the original variables: 

xi{t) = -x 2 {t) = -\yi{t) 
ui = -\yx- \V2 ~ f 2/4 
u 2 = \yi - \y2 , 

with xi(io) = —32/1(^0) a given intensity at time to- 

It is clear that the current intensities, which have almost surely continuous paths, are much 
more regular than the voltage drops, which are only random distributions. 
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